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ABSTRACT 

We propose that supermassive stars may form in quasar accretion disks, and we discuss possible 
observational consequences. The structure and stability of very massive stars are reviewed. Because 
of high accretion rates, quasar disks are massive and the fringes of their optically luminous parts are 
' prone to fragmentation. Starting from a few hundred solar masses, a dominant fragment will grow to the 

I isolation mass, which is a significant fraction of the disk mass, more quickly than the fragment contracts 

. onto the stellar main sequence. A gap will form in the disk and the star will migrate inward on the 

accretion timescale, which is comparable to the star's main sequence lifetime. By interrupting the gas 
, supply to the inner disk, the gap may temporarily dim and redden the quasar. The final stages of stellar 

D ' migration will be a strong source of low-frequency gravitational waves. 

. Subject headings: accretion disks — gravitation — quasars: general 

o : 

^ ■ 1. INTRODUCTION 

CNJ ' Accretion disks massive enough to fuel bright quasars are expected to be self-gravitating beyond a few hundred to a few 

^ ' thousand Schwarzschild radii (i?s), {e.g. Shlosman & Begelman 1987) Even under extreme assumptions, it is difficult to 

, see how such disks could extend beyond ^ 0.1 — Ipc without fragmenting completely into gravitationally bound objects 

' (Goodman 2003; Sirko & Goodman 2003). Self-gravity is less problematic beyond ~ 10 pc because the stellar bulge 

' dominates the rotation curve, and because the disk eventually becomes optically thin. 

: Quasar accretion disks are not resolved and there are few direct constraints on their structure. The principal reason for 

' believing in them at all is that one does not know of any other plausible mode of accretion that converts mass to radiation 
with high efficiency. The overall size of the accretion flow is however constrained by the need for sufficient surface area to 

■~~>~ radiate the observed continuum. At wavelength A, this area is 

^ ■ . -1 . . 3 
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assuming incoherent emission by gas at temperature Tg. Thus the optically luminous part of the disk should extend to 



^ ■ at least ~ lO'^iJgj where self-gravity becomes important. 

I We shall assume that a conventional thin disk does exist and extends to 10'^i?s with an accretion rate sufficient to 
J> , support a typical bright QSO of mass > 10® Mq and luminosity > 10^^ ergs~^. We shall argue that a likely consequence 
of the incipient self-gravity at the outer edge of such a disk is formation of supermassive stars. 

The term "supermassive" is truly justified here. While fragments may start with only IO^Mq, they seem likely to 
' grow to a significant fraction of the disk mass, perhaps ~ lO^M©. Yet it is quite possible that even such extreme objects 
would escape notice until appropriate gravitational- wave detectors are built (§5). 

Supermassive stars of even larger masses (~ 10®) were once proposed as models for quasars themselves (Hoyle & 
Fowler 1963; Zel'dovich & Novikov 1971). This notion has long been abandoned on grounds of stability (§2), and more 
conclusively because statistical arguments demonstrate that quasars convert mass to energy much more efficiently than 
nuclear fusion (Soltan 1982; Chokshi & Turner 1992; Yu & Tremaine 2002). It remains possible that supermassive stars 
may be seeds or precursors to quasars (Rees 1984). 

In this paper, we accept the conventional view that the energy source for quasars and AGN is accretion onto massive 
black holes. But the disks required by this interpretation of the facts are so massive and dynamically cool that they are 
likely to form substructures. We expect this to occur via the standard local dynamical instability of self-gravitating disks, 
i.e. the gaseous version of the Toomre instability (Toomre 1964), although growth of very large masses may perhaps be 
seeded by other means, for example capture of stars (§4.3). 

In steady accretion, gravitational instability becomes progressively more severe towards larger radii {e.g. Goodman 2003, 
henceforth Paper I). Thus one might chose to consider a strongly unstable region where the disk is likely to fragment 
completely: in other words, where most of the mass would be found in clumps rather than a smooth layer. Others have 
investigated this regime (Shlosman & Begelman 1989; Shlosman et al. 1990; Kumar 1999). However, we have chosen 
to focus on a transitional region where self- gravity is mild and bound structures, where they exist, are likely to be well 
separated. There are several reasons for our choice. For one thing, we doubt whether very strongly self-gravitating 
material can form an accretion flow at all, except as stars (Paper I). For another, there is little direct evidence that steady 
disk accretion extends beyond ^ 10'^i?S: where self-gravity becomes strong. If it does, then a marginally self-gravitating 
region surely exists by continuity. Finally, marginal self-gravity offers the methodological advantage that one can study 
it as a perturbation to non-self-gravitating, thin-disk accretion, which is somewhat understood. 
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The plan of the paper is as follows. In §2, we review the structure of supcrmassive stars and what is known about their 
stability. §3 describes the physical conditions of a QSO accretion disk in the region where it is likely to be marginally 
self-gravitating, assuming an a disk in steady state (Shakura & Sunyaev 1973). Scalings with black hole mass, accretion 
rate, and viscosity parameter are given. §4 builds upon these standard results to discuss the formation, growth, and fate 
of bound fragments ("stars") that form in marginally self-gravitating regions, with particular emphasis on their maximum 
mass. §5 sums up. 

2. VERY MASSIVE STARS 

This section collects some elementary but important structural formulae for very massive stars; somewhat more accurate 
forms were given by Bond et al. (1984, henceforth BAG). We also summarize what is known about the stability of such 
objects, both before and after they reach the main sequence. 

Very massive stars tend to be radiation-pressure dominated. To the extent that they are chemically homogeneous and 
the opacity (k) is constant, in hydrostatic and radiative (but not necessarily nuclear) equilibrium, they are approximated 
by Eddington models: 

Pgas = P*P + constant, 
L, = (l-/3,)LEdd,* = (l-/3*)^^^^^, (1) 

K 

in which pgas and prad = P — Pgas are the partial pressures of gas and radiation, and the asterisk is used to distinguish 
parameters of the star from those of the surrounding disk. The dimcnsionlcss parameter /3* < 1 is nearly constant within 
a given star. The additive constant on the first line above is unimportant throughout most of the interior, so that we may 
regard /3* approximately as Pgas/P- 

Except near the surface where the additive constant must be taken into account, and in the core where luminosity may 
not be proportional to mass, it follows that the stars are n — 3 polytropes, 

-,1/3 
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where is the molecular weight relative to hydrogen, so that fj, w 0.62 for a fully ionized gas of solar metallicity. Prom 
standard results for such polytropes, 

''*~47^f^V- (3) 



It is important that this relation does not depend upon the stellar radius (i?*), so that it holds even if the star has not 
reached the main sequence. Thus for extreme masses M ^ lOOM©, one has /3* oc M~^/^. The binding energy 

is small because the gravitational and internal energy nearly cancel, and the Kelvin-Helmholtz time is correspondingly 
reduced: 

- ISTTCi?. 1-/3/ ^ > 



The entropy per unit mass is 

S- 
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+ constant, (6) 



so that the effective adiabatic index for perturbations is 



_ 32 - 24/3. -3/3; ^4 1 2^ 

On the main sequence, the central temperature is nearly constant. Since GM^/R^ a Pc T^/{1 — P*), one has 
iJ* oc My^. More accurately (BAG) 

D n ^0.39 / Y \ 0-05 / M \ ° / „ \ 

ii*,ms _ ir,\^ ~ P*> I -^CN \ , i c / -'^* I ( 1^* \ ^ ir,3i 



where Xcn is the mass fraction in GNO elements. The Kelvin-Helmholtz time on the main sequence is therefore 

tKH,ms« 3300 ( ^ ( ^V/3-""'^-'(l-/3,)-"«Vr. (9) 



^0.4gcm~2y \^ ^ 

The dependence on M*, which enters through /3. and eq. (3), is so weak that we shall generally ignore it. 
The main-sequence lifetime is also nearly constant (BAG), 

tms w (2 - 3) X 10'' yr ~ — , (10) 

-t^Edd,* 

unless cut short by catastrophic instabilities. 
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2.1. Stability of massive stars 

Instabilities can be categorized according to if they occur before, during, or after the main-sequence phase. 

The most famous pro-main-sequcncc instabihty is a rclativistic one. Since the adiabatic index (eq. 7) is close to 4/3, 
changes in internal and gravitational energy are almost equal and opposite if the star is homologously and adiabatically 
compressed. Small corrections to the energy can have a large effect on stability. One such correction is general relativity, 
which increases (in absolute magnitude) the potential energy by a fraction ^ (GAI^/c^R^). A very crude estimate of 
the minimum unstable mass can be obtained by evaluating GMsf/R^f on the main sequence and equating this to j3<t/&- 
Using equations (8) and (3), this yields M^^qr ~ 6 x 10"" M0. More detailed calculations predict a rather smaller mass 
(Chandrasekhar 1964), M^^^gr ~ 5 x IQ'^Mq, for nonrotating stars. 

Rotation lends stability because its energy scales with pgas rather than Prad under homologous contractions, that is, 

Ey-ot cx R~'^ oc pV^ ■ The maximum rotational energy of a uniformly rotating n = 3 polytropc is small, i^rot/l-Ej < 0.007, at 
which value centrifugal and gravitational forces balance at the equator (Zel'dovich & Novikov 1971). General-relativistic 
effects are also small, however, so rotation is important, especially since stars forming out of a disk will probably rotate 
strongly. We therefore take the minimum mass for general-relativistic instability to be the value appropriate to maximal 
uniform rotation (Baumgarte & Shapiro 1999) 

M,,GR « lO^M©. (11) 

Stars with M* < reach the main sequence on a timescale no longer than the Kelvin time (cq. 9). Thermonuclear 

reactions commence on the CNO cycle, since quasar accretion disks are at least as metal rich as the Sun (Dietrich et al. 
2003, and references therein). At > IOOM0, the sensitivity of the thermonuclear reaction rate to central temperature 
drives a linear ovcrstability to pulsations (Eddington 1926; Lcdoux 1941; Schwarzschild & Harm 1959). Non-linear 
calculations suggest that this instability drives mass loss (Appcnzellcr 1970b, a; Papaloizou 1973). The rate of mass loss 
is uncertain, but is probably no greater than m* = 5 x IO^'^^Mq yr^^ for a 13OM0 star and m* = 5 x 1O~^M0 yr^^ for 
a 27OM0 star (Appenzeller 1970b), and could be as little as ~ 1O~^M0 yr~^ for similar masses (Papaloizou 1973). The 
scaling of to* to much larger stellar masses is also uncertain. 

The above mass- loss rates may be compared to those resulting from radiation pressure driven winds. Consider a I2OM0 
zero age mam sequence (ZAMS) star with T^s = 5.3 x 10'' K and luminosity L* = 1.8 x 10'^ Lq (Schaller et al. 1992). 
The analytic relation of Vink et al. (2000) for radiation pressure driven winds including multiple scattering, predicts a 
mass-loss rate of 1.65 x 1O~^M0 yr~^. This mass-loss rate scales approximately as L*'^ and to~^'^, with an additional 
weaker increase resulting from hotter surface temperatures. It can be considered as a lower limit to the actual mass-loss 
rates from massive stars, since these typically increase during the later periods of stellar evolution, particularly during 
the Wolf-Rayet phase. For example, at the above rate the stellar mass would decrease to about 7OM0 in the lifetime of 
~ 3 Myr. However, more detailed evolutionary models including mass loss rates calculated consistently for each stage of 
the evolution give a final stellar mass of only ^ 8M0 (Schaller et al. 1992); half of the total mass loss occurs in the last 
0.5 Myr of evolution. It should be noted that these results are for non-rotating stars, evolving with no additional mass 
accretion. 

If the star manages to survive the main-sequence, then in the post-main-sequence phase an instability arises by the 
creation of e+e~ pairs when the central temperature ^ 10^ K (Zel'dovich & Novikov 1971). This is believed to result in 
an explosion if M* < 2OOM0 (BAG) but complete collapse to a black hole if M* > 300Mq (Fryer et al. 2001). 

3. INITIAL CONDITIONS FOR MASSIVE STAR FORMATION 

Several physical properties differentiate the star-forming environment of a quasar disk from that of a typical region of 
Galactic massive star formation (Plume et al. 1997; McKce & Tan 2003). Firstly the densities and pressures are many 
orders of magnitude greater. Secondly the environment is subject to strong shear, due to the Keplerian orbits around 
the massive black hole. The dominant source of heating of the gas, at least in the inner regions of the disk, is due to its 
accretion in the QSO's potential. At the typical temperatures in the QSO disk we expect dust grains and molecules to 
have been destroyed, and finally, radiation pressure can lend significant support to the gas. 

3.1. Physical conditions in the disk 

From Paper I, we have the following properties for geometrically thin, optically thick, Keplerian disks that are heated 
by viscous dissipation only. The radial dependence of the surface density, S, and the midplane temperature, T, are given 
by: 

/ K^miH \ ^2/5j^3/5 (12) 

« 5.27 X 10^ ( , \ M^'/'rf/'° K, 

o4/5 / 4 4 \ -'-/^ 

S = |_ / t^J^ \ (a/?''-l)-4/5K-V5M3/5fi2/5 (13) 

« 2.56 X 10^ {ao.sp'-'r^/H'J'^ol^'^-'^'f^^^'M'^/'r-'/' g cm'^ , 
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The meaning of the various symbols is as follows: fl — (GM/r^)^^^^ , k is the opacity relative to the electron-scattering 
value (0.4 cm^ 1^ is the mean particle mass in units of the Hydrogen atom mass, niH', O.Sap.s is the Shakura-Sunyacv 
viscosity parameter [scaled to the largest likely value in a selfgravitating, but not fragmented, thin disk (Gammie 2001)]; 
(3 = Pgas/pi 0-l£o.i is the radiative efficiency of the disk {— L/Mc?)\ is the disk luminosity relative to the Eddington 
hmit; M = WM^Mq is the black hole mass; and i?s = 2GM/c^ w 10 ^Mgpc is the Schwarzschild radius, with the 
orbital radius in units of lO^iJg- Also b takes a value of either or 1, depending on whether the viscosity is proportional 
to the total or gas pressure, respectively. 

By including the opacity as an explicit parameter, we avoid its sometimes complicated functional dependence on density 
and temperature, thereby obtaining simple algebraic formulae for the disk properties. The numerical disk models of Sirko 
& Goodman (2003, henceforth Paper II), which took realistic opacities, show that k ~ 1 gcm""^ throughout the region of 
the disk of interest to us, namely r < 10^i?s- (Indeed, k ~ 1 out to r 10'' -Rs if nonviscous auxiliary heating maintains 
Q » 1 to that radius.) Furthermore, as will be seen, many quantities of interest depend rather weakly on k. 

If viscosity scales with gas pressure {b = 1) then eqs. (12) & (13) do not depend on /3, which in any case is a known 
function of density and temperature: 

P _ Pgaa _ 3cA:b _P_ _'^( 



using eqs. (12) & (13), 

^(l/2) + (b-l)/10 



(237r4)i/5a-i/inc(A;B/m)'/'^"'/'°«"^/'°i^"^/'°M-4/5 



1-/3 

« 0.311 a-J/'°/c-^/iV-'/'(£o.i/«K)'/'M-^/^° (14) 



Using eq. (12) to eliminate the isothermal sound speed, Cg = \Jpl p = kT/{iimH) from Toomre (1964)'s gravitational 
stability parameter yields 

Q = 3(47r)-3/5a7/i0/3(''-''^/'° ( -^V^' ^^''^'^G-Ik^/iom-^/^j^s/io 

\HmHj 

« 1.0 ao^' V^'''-''^/'°eo^'/E '/'«'/'V-'^/'M8-^^/^°r3^^/^° . (15) 

Consider quasar accretion disks at the inner radius of gravitational instability where Q = 1: 

r,.g. « 915 aJ^/''«2/V-«/^(eo.i/;E)'/''M-'*^/''/3(i4''-24)/27i?s. (16) 

This equation differs slightly from cq. (10) of Goodman (2003) because some approximations assuming /3 <C 1 were used 
there. 

Substituting for r from eq. (16) into cq. (14) yields 

ol -6/3 

1^)3/4 = 0.389 a;/3^ft-VV-H^o.i//.)^/^M-^/« . (17) 

Thus gas and radiation pressure are comparable at rg.g.. Equations (17) & (16) imply rg.g. is not very sensitive to the 
choice of viscosity law, e.g. r^.g ~ 2700i?s for b ^ vs. 1700i?s for 6 = 1 at the fiducial values of the other parameters. 
On the other hand, radiation pressure is at least moderately dominant at all r < Tg.g., so that we may replace 1 — /3 with 
unity in the quartic (14); this makes /3 a simple power law in all other parameters. To further simplify the discussion, we 
shall also take the viscosity to be proportional to gas pressure except where stated otherwise. 

With these assumptions (/3 <C 1, b = 1), the properties of the disk at r < rg.g. are as follows. The midplane density is 

p=^^ 4.31 X 10-i°ao.f «-«/V/'(eo.i/?B)^/^M-^/^3^/^gcm-3 (18) 

(corresponding to number densities of hydrogen nuclei ^ 10^^ cm~'^), and the total pressure is 

p/fcB « 1.41 X 10^0 a;;/|K4/5^^-4/5(;^/^^^)8/5^^4/5^-18/5 ^^-3 ^ (^g) 

The midplane temperature is given by eq. (12) with 6=1 (note that at this temperature dust and molecules will be 
destroyed), which corresponds to an isothermal sound speed, Cg = \ppfp^ of 

Cs w 67.2 /tV5/^£-ir-3"^/^° kms"^ . (20) 

The disk scale height is 

/i=g«10K/K%iiis, (21) 
the surface density is given by eq. (13), and the total disk mass inside r is 

Mdisk(< r) « 5.03 X 10^ or^^^k-^l''p^l\l^le^.^fl^Mfl\y^MQ . (22) 

The orbital timescale is 

torb« 8.79 Mgr^/^r, (23) 
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whereas the accretion time is 



t^,,^ = M^^iE^A « 2.28 X 10^ a^',^'k-'/'n^/Heo.i/Lr/'Ml/V,/' yr , (24) 



Wcirc w eriOrg ^/^kms-^ (25) 



and the circular velocity is 

„, . ^ «7i n y ^ 

Finally, we may eliminate (3 from eq. (16) to express the radius of marginal self-gravity itself as 

rs.g. « 1550(ao.3//i)'/'(/E/£o.i)'/''«'/'M8"'/' Rs- (26) 
It is worth emphasizing once again that eqs. (18)-(26) assume u oc Pgas and /3 <C 1 at r < rg.g.; in view of eq. (17), the 
latter condition is expected to be satisfied at high accretion rates, M >! Mq yr~^. 

3.2. Fragmentation 

We assume that the initial sizes and masses of protostars forming in the disk are determined by Toomre (1964)'s 
dynamical instability. The shear stabilizes long wavelength perturbations from collapse, while the pressure stabilizes 
short wavelength modes. The most unstable mode has radial wave number fc^.u. = TrGS/c^ = {Qh)~^ so that 

27r 

Am.u.(^s.g.) = T ~ 63K(/E/eo.i)^s- (27) 

"-m.u. 

The numerical expression assumes /3s. g. <C 1 so that h is given by eq. (21), which is probably a good approximation for 
Mg > 1 and Ze ~ 1- The instability is axisymmetric, but fragments having comparable azimuthal and radial dimensions 
probably result. The corresponding mass at r = rg.g. is then of order 

M^,^, = EA^.,. - 300 a^\fi'/^ii{l^Ms/eo.if'''MQ. (28) 

A second way to estimate the fragment mass is to note that (3 is approximately conserved under an adiabatic contraction, 
since as eq. (6) shows, the entropy depends only logarithmically on density at fixed /3. Therefore one may set P{rs.g ) 

in the relation oq. (3) for Eddington stellar models. This yields an initial stellar mass about four times larger than eq. 
(28) but with the same scalings. The star is likely to gain a great deal more mass by accretion, as we shall explain next. 



4. ULTIMATE MASSES 

In a strongly unstable case where Q < 1, the disk might break up completely into closely spaced fragments. The 
subsequent evolution would then be dominated by collisional agglomeration. Let us consider instead a less unstable 
situation (Q 1) so that the initial fragments arc well separated and contain only a small fraction of the disk mass. Then 
a protostar gains mass mainly by accretion from the ambient gaseous disk. 

We first consider whether the fragment (with mass given by eq. 28) will itself fragment into smaller masses. The 
usual condition for this to happen is that the cooling time be less than the the dynamical time. The initial density 
(p,) and optical depth (r* ~ Kp,i?,) of the fragment will be comparable to those of the ambient disk, and its radius 
R* ~ Am.u. ^ 2/i, the disk half-thickness. The initial dynamical time is (Gp*)^^/^ ^ (Gpdisk)^^^^ = and the cooling 
time of the fragment is of order the local thermal time of the disk, t^h = a~^^ (Pringle 1981). Therefore, provided aQ < 1 

then icooi > idyn- Furthermore, if the fragment contracts at constant mass, then tdyn oc R*^^, whereas icooi oc R^^ if k 
is constant, so that tcooi/idyn oc R7^^'^, making further fragmentation increasingly unlikely as the radius shrinks. Even if 
subfragmentation were to begin, rapid accretion from the surrounding disk (see below) would probably drive the fragments 
together again by absorbing energy and angular momentum from the relative motion of the fragments. In short, the mass 
of the protostar is more likely to grow by accretion than to decrease by fragmentation. Nevertheless, it would be useful 
to confirm this by numerical simulations in two or three dimensions with allowance for cooling. 



4.1. Accretion and the isolation m,ass 

The accretion rate will be estimated by analogy with growth of terrestrial planets in a planetesimal disk (Lissauer 1987). 
An important lengthscale in the latter problem is the Hill radius Rh = (Afp/3M,)^/^r, which is essentially the size of 
the Roche lobe surrounding the planet (mass Afp) at an orbital distance r from the central star (mass Af,). At distances 
<C Rh from the planet, motions of planetesimals are dominated by the gravitational field of the planet rather than the 
star. In a quasar disk, the role of the planet is played by the protostar (mass M* 1O^~^M0), and that of the Sun, by 
the central black hole (mass M > IO^Mq). Thus 

«».(§)"%. (29) 

Tidal torques exerted by the growing protostar will eventually open a gap in the disk, shutting off or at least reducing 
the rate of accretion. These torques, however, act only on gas that has encountered the protostar at least once. Material 
on an orbit whose semimajor axis differs from that of the protostar by ~ Rh stands a good chance of being accreted on its 
first passage. Let us therefore assume that the protostar accretes the entire annulus |r — r*| < /hRh, where fn ^ G(l). 
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(Lissauer infers fn 3 — 4 from simulations of planctesimal growth.) The mass of this annulus is larger than the initial 
mass (eq. 28) of the protostar by a large factor ~ r/ft. if Q ~ 1. Therefore, the initial width of the gap that is cleared is 
proportional to Hill radius of the mass accreted. This condition defines the isolation mass (Lissauer 1987), 

= ^'Y^'^lfJ^' « 0.96 X 10^ /^/^ao.3/'/?''^^-'^/'^-'/^V^/^(«./^o.i)^/^°M3"/Vf /^°Mo. (30) 

Numerically, if Jh ^ 1, this is comparable to the disk mass (eq. 22), but formally Miso oc Mdisk x (M^isk/My^^ . The 
dependence of Mjgo with r and M is shown in Figures 1 and 2. The isolation mass is startlingly large and one wants to 
examine very critically the assumptions that have led to it. 

The dynamics of a pressure-supported gas are somewhat different from those of a nearly collisionless swarm of plan- 
etesimals. Rafikov (2003) has recently demonstrated that the growth rate of a planetary embryo slows significantly at 
a mass <^ Miso because gravitational encounters with the embryo tend to "heat" the epicyclic velocity dispersion of the 
planetesimals and thereby decrease their rate of collisions with the embryo. In this respect a quasar disk is more favorable 
for rapid growth because the gas can cool radiatively. In fact, the cooling time is ~ {afl)~^, which is probably short 




30.01 0.1 1 10 102 103 0.01 0.1 1 10 102 103 



r3 = r/103Rg 



rg = r/l03R, 



Fig. 1. — Isolation mass, Miso (sec cq. 30), versus radial location (in Schwarzschild radii) in the disk (heavy lines) with 6 = 0, 1 (bottom, top 
panels) for lO'', 10*. lO'^Mg eentral black holes with k = 1, a = 0.3, /i = 0.6, Ie = eo.i = li and fn = 1 (left panels). In each panel the lower 
thin line shows the locus defined by r = r-s.g. (eq. 16), with each point on the line representing a different central black hole mass. Similarly, 
the upper thin line shows the condition when Migo(r) = M^ig^.{r) (here we define M(jig]j(r) = (10/7)7rS^ as a measure of the maximum mass 
available for star formation from the disk). The dotted line traces the conditions where the midplane temperature is 10* K — in the cooler 
region we may expect significant departures from the fiducial opacity. The right panels show the effects on Miso(r) for the Aft = 1 case of 
increasing k by a factor of 10 (heavy dashed line) and reducing a by a factor of 10 (heavy long dashed line). The effects on the self-gravity 

3/2 

radius condition and the disk mass condition are also shown with the same line types. Note that Mjao oc , which introduces additional 
uncertainty. We expect supermassive star formation to the isolation mass at r ~ ''s.g.i and from these figures we can see that this is not 
prevented by a lack of material from the disk for the range of parameters we have considered. 



T \ I I I I I r 




M (Mq) 

Fig. 2. — Isolation mass at r = rs.g (top panel) and lO^Rs (bottom) versus central black hole mass, assuming a = 0.3, « = 1, = £o.i 
and fn = 1- The cases with 6 = 0, 1 are shown, as labeled. The dashed lines show the disk mass inside r. 
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compared to the time between successive encounters with the protostar, ^ 27r(f2///i?///r)^^. Also, the accreting protostar 
probably nearly fills its Hill sphere during the growth phase (see below) so that it presents a relatively large cross section 
for "collisions" compared to a solid planetary embryo whose density is determined by chemical bonds. 

On the other hand, whereas impacts of planetesimals on planetary embryos are perhaps fully inelastic, it is not imme- 
diately obvious that the accreting protostar can accept gas at the rate that it arrives via differential rotation. To address 
this question, we compare the rate at which mass enters the Hill sphere, the rate at which this material sheds its "spin" 
angular momentum, and the rate at which it sheds energy. Since masses much larger than (28) are concerned, we estimate 
these rates in the limit Rh ^ H. 

We begin by estimating the rate of growth of the Hill sphere. Material on an orbit at radial separation Ar from the 
center of the Hill sphere approaches it at azimuthal speed 

Vrei = r—Ar w --ilAr. 
dr 2 

The surface density on such a streamline is S, the background surface density of the disk, if the material is approaching 
the star for the first time, that is to say, 

27rr 47r r ^ 1 ,„ , 

where t is the time since accretion begins. We have set Ar = /hRh because if the inequality is satisfied at this separation, 
then it is satisfied at smaller separations, and more distant material is not likely to be much perturbed. Hence mass enters 
the Hill sphere at the rate 

r+fnRH 

Vrel\dr ^ -Y.\1}hRh- 



dt J ' ' 2 

r—fHRn 

Since Rr = r{MH/SMy/^, this can be rewritten as 

—Rh = -^TTT = constant. (32) 
dt 6M ^ ' 

(We use Mh rather than M* to denote the mass within the Hill sphere so as not to prejudge the question as to the 

timescale on which this mass concentrates into a quasi-spherical star rather than a centrifugally supported disk: see 

below.) Thus Rh increases linearly with time until the two sides of the inequality (31) become equal, at which point 

accretion presumably slows or stops. This defines the asymptote of Rh and hence the isolation mass: 

'27rSr2\^/^ 



i?H(M.o)=^^j r, 
_ (27r/«Sr2)3/2 



9MV2 

The mass is smaller than our previous estimate (30) by a factor 6^"^/^, which is an indication of the crudcness of both 
arguments. The approximate time at which this mass is achieved is found by interpreting eq. (31) as an equality with Rh 
as given above: 



27r / 2M ^ 



«210/^-^/^«^/3^«i/i"M-2/5(£o.i/ZB)-^/^"Ms2/Vf yr. (33) 

Since M/Y,r^ ~ r/H at Q = 1, tiso is longer than the local orbital period by ~ {r/Hy^^. 

Having estimated the timescale on which the isolation mass accumulates within its Hill sphere, we now consider the 
timescale on which it contracts to radii <^ Rh- Unless the contraction time is smaller than the accumulation time, it is 
questionable whether the Hill sphere can accept the full mass Miso in the first place. On the other hand, a contraction 
time much shorter than tjso would suggest fragmentation. 

In order to contract, the mass Mh must shed both energy and the angular momentum associated with rotation around 
its own axis i.e. spin. We will now show that the timescale for loss of angular momentum is likely to be shorter than that 
for loss of energy. Suppose the contrary: then the material within the Hill sphere settles into a centrifugally supported 
disk, which will be called the "stellar disk" to_distinguish it from the larger disk that orbits the black hole. Let the 
mass-weighted average thickness of this disk be h, presumed to be ^ Rh but not necessarily equal to H, the thickness of 
the ambient black-hole disk; let its median (half-mass) radius be R; and let its median angular velocity around its own 
axis be u. Then the viscous time in the stellar disk is {au))~^{R/hy . Suppose provisionally that this is long compared 
to tiso, as it would be if i? ^ Rh and h < H. Then R cx Rh and <D oc 17 as Rh expands, since the specific angular 
momentum of the material falling into the Hill sphere is ^ ilR^. [In fact, the stellar disk should be no larger than the 
streamline-crossing radius, which is i?sc ~ O.ARh at these extreme mass ratios Mh ■ M_< 10^^ (Paczynski 1977). At and 
beyond i?scj angular momentum is efficiently removed by the tidal potential.] Thus h oc Cg/fi, where Cg is the median 
sound speed in the stellar disk. The cooling rate of the disk by vertical radiative diffusion is therefore 
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tKn.H « 53.0ao.f «i^/2V-^/^aB/£o.i)^/^"M-^/Vf (^\ yr, (36) 



where /i is a positive dimensionless constant, c the speed of Ught, k the opacity, and S the median surface density. We 
have assumed that the stehar disk is radiation-pressure dominated, hke the ambient black-hole disk, but if not, then the 
cooling time would be longer by a factor (1 — 

New material joining the disk does so through a shock that endows it with internal energy a GMh/R per unit mass. 
Thus the median internal energy per unit area of the circumstellar disk evolves as 

in which /2 is another positive constant. From the discussion above, R oc Rh oc t, Mh oc t^, and therefore Mh cx; and 
E oc Mh/R% oc t. So if h were constant in time, then would also be constant, and the three terms in eq. (35) would 
scale with time at different rates: the first two would be constant, but the third oc t^. Thus the equation could not hold 
for all t. The solution is to let Cs oc t so that the first and third terms balance, in which case h oc Rr- 

In short, a ccntrifugally supported stellar disk would not remain thin but would evolve toward a quasi-spherical, 
pressure-supported configuration as long as it is rapidly accreting from the ambient black-hole disk. Once h ~ Rh, the 
viscous time becomes ~ (aJl)^^ at most. In fact, if the "star" fills its Hill sphere, then the angular momentum of its 
outer parts is not conserved, while its inner parts may transfer angular momentum outward by global bar instabilities on 
their own dynamical timescale, which is < Q.~^ . 

At this point, it appears that energy is more problematic than angular momentum. The estimate (34) suggests that 
the cooling time of a quasispherical configuration should scale as t^, whereas Rr/Rh oc t. This would seem to imply that 
cooling cannot keep pace with accretion, and therefore that the star should fill or overflow its Hill sphere. But the matter 
is more subtle than this because of the peculiarities of F w 4/3 polytropes. 

For a completely nonrotating star of fixed mass, the contraction timescale is the Kelvin-Helmholtz time obtained from 
eqs. (5) and (3); putting i?* = Rh, 

1/6 

\M~ 

which is somewhat shorter than the accumulation time tiso [cq. (33)]. (Henceforth we write M» rather than Mh for the 

— 1/2 

mass within the Hill sphere.) The critical difference between the two timescalcs is the factor oc M« in eqs. (4) & 
(5), which reflects the nearly zero binding energy of a spherical high- mass star, so that only a small fraction of the internal 
energy need be radiated in order to cause the star to contract substantially. This factor does not occur in the vertical 
contraction of a disk, of course. 

Furthermore, the star could contract even without radiative losses because it accretes material with negative total 
energy. It can be shown that when Rh ^ H, so that fluid streamlines can be approximated by ballistic test-particle 
trajectories, material accretes onto the Hill sphere with a negative Jacobi constant per unit mass, 

where is the center of the star, Q. its orbital angular velocity, and v is the test-particle velocity relative to the 
star in a frame rotating at $7. (The Jacobi constant of the accreting streamlines cannot be arbitrarily negative, since 
Hj > — (9/2) (ni?i/)^ in order to cross the inner or outer Lagrange point.) Very close to the star, at |r — r*| -C Rh, Hj 
reduces to the usual kinetic-plus-potential energy in the potential of an isolated mass M^,. Thus, the putatively spherical 
star accretes material that has negative total energy even without radiative losses. This would cause contraction on a 
timescale even shorter than (36) if rotation could be ignored altogether. However, we expect that the star will in fact 
rotate with a median angular velocity scaling in proportion to Vt as it gains mass. Consequently, the star's rotational 
energy cannot be entirely neglected unless it is smaller than the binding energy given by eq. (4), which is marginally 
unlikely since ~ 10~^ at Mjso. 

To recap, the growing star cannot cool fast enough to become a disk but will be more nearly spherical and therefore 
able to shed angular momentum quickly; on the other hand, the residual rotational energy, though small compared to its 
gravitational energy, seems likely to stabilize it against rapid contraction. The interplay between energy loss, angular- 
momentum loss, and radiation-pressure dominance (so that F w 4/3) is so intricate that it will be difficult to decide 
whether the star overfiows or detaches from its Hill sphere without elaborate calculations that are beyond the scope of 
this paper. 

Again, it would be useful to have recourse to numerical simulation. We are not aware of any calculations for disks 
supported largely by radiation pressure, but Bate et al. (2003) and Lufkin et al. (2003) have recently simulated three- 
dimensional accretion from a protostellar disk onto embedded planets ranging in mass from that of earth (Mg) to that 
of Jupiter (Mj). Both studies place their planets at rj = 5.2 AU and follow the evolution for ~ 10^ orbits in a disk with 
h/r = 0.05. Bate et al. use an eulerian finite-difference method, whereas Lufkin et al. use smooth particle hydrodynamics. 
Neither includes cooling, but this is not an impediment to accretion because the former study removes removes the gas 
at each time step from the grid cells closest to the planet (well within the Hill sphere), while the latter uses an isothermal 
equation of state and includes the sclfgravity of the gas, which allows indefinite compression of the planet's atmosphere. 
The accretion rates measured in these simulations are indeed comparable, up to factors ^ 2, to what is predicted by eq. (33) 

with /h = 1, even though the most massive planets may exceed the isolation mass (eq. 30), Mi,o = (0.089, 0.25) /^/'Mj, 
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corresponding to the surface densities S(rj) « (0.23, 0.46) Mj/rj in the two studies. (Bate et al. impose a partially 
cleared gap around the planet in their initial conditions; wo have quoted S as it would be if there were no gap.) 

Unfortunately for the purposes of this comparison, the Hill radii of the most massive planets in these simulations is 
scarcely larger than the disk thickness, so that the accretion is not very far into the two-dimensional regime contemplated 
by equation (33). To express this another way, the planetary isolation mass quoted above is comparable to the mass 
defined by RniMh) = h, viz. Mh = 3{h/r)^MQ « 0.38Mj. The ratio M^/Mso scales as {Qh/rf/"^, which is of order 
unity for gaseous protoplanetary disks but very small at r ~ 10^i?s in quasar disks because both h/r and Q are smaller 
there. Satellites ^ Mh probably cannot open a gap regardless of the viscosity of the disk (Lin & Papaloizou 1993); they 
will undergo Bondi accretion provided that the accumulated gas cools rapidly enough. Thus, the applicability of the 
concept of isolation mass in quasar disks remains to be tested. 

Once the isolation mass is achieved, or possibly even earlier, the star will detach from its Hill sphere and contract until 
it reaches the main sequence after a time of order the Kelvin time (cq. 9). It is important to note that, at least for our 
nominal parameters, the tim,e required to reach the main sequence is longer than the time to accrete the isolation rn,ass: 
therefore we do not expect nuclear burning to impede accretion. The formation time is also short compared to the viscous 
evolution time of the disk (eq. 24). 

4.2. Orbital migration and tidal disruption 

When fusion begins, the star may disintegrate via pulsational instability (§2.1). It would certainly be interesting to 

understand this stage: will the star disrupt completely or leave behind a remnant? How much will the disk be enriched? 
Since the answers to these questions are not available at present, we assume provisionally that the star survives the 
instability with much of its mass intact. If so, then it can be expected to undergo Type II radial migration: that is, while 
maintaining an annular gap in the disk by tidal torques, the star will drift inwards on the viscous timescale (eq. 24) by 
analogy with jovian planets in protostellar disks (Ward 1997). For our nominal parameters, the latter timescale is shorter 
than the main-sequence lifetime (cq. 10). Thus the migration may be completed before the star leaves the main sequence. 
The star will be tidally disrupted at an orbital radius where Rh equals the main-sequence radius (eq. 8), 

/ \ 0.95 / ,^ \ 0.14 

which is comparable to the radius of the marginally stable orbit for a Schwarzschild black hole with M = IO^Mq. 

4.3. Capture 

The concept of isolation mass is independent of the initial Q of the disk, although it is true that Mi^o increases with the 
overall mass of the disk [eq. (30)]. Thus mechanisms other than gravitational instability may provide the initial seed from 
which the isolation mass develops by accretion. In the primordial solar nebula, for example, where most likely Q ~ 10^ at 
the present position of Jupiter, that planet is believed to have been initiated by coUisional agglomeration of solid bodies. 
We do not expect solids in the non-self-gravitating parts of quasar disks, but there may be other ways to form seeds that 
might grow to the isolation mass. One such process is capture of pre-existing stars in collisions with the disk. 

The capture process has been examined at length by Syer et al. (1991) and Artymowicz et al. (1993), and we have only 
a few remarks to add. In the following discussion, all disk and black-hole parameters are assumed equal to their nominal 
values in §3 except the orbital radius, for which we assume rs < 1 so that self-gravity is slight and radiation pressure 
dominant. 

The first remark is that in these regions, stars may perhaps more often be destroyed by their encounters with the disk 

than captured. The energy that must be dissipated to circularize a typical stellar orbit exceeds the binding energy of a 
solar-type star by a factor ~ lO'^rg ^. Most of the dissipation should occur in a bow shock driven into the tenuous disk 

gas rather than within the star. From §3.1, the peak postshock pressure is ~ pv'^ ~ lO^r^^^^ dyne cm~^; at just ~ 10~^ 

of the star's central pressure, this is only a gentle squeeze. However, except near the stagnation point, the postshock flow 

passes the star at ~ lO^r^^^'^ kms~^ and at a temperature between lO^r^'^^^ K and ~ lO^r^^^^ K (depending whether the 
radiation field comes to equilibrium with the plasma; the mean free path for creation of soft photons by bremsstrahlung 

is comparable to the standoff distance of the shock), so that it may erode the stellar surface, which is not rigid but simply 

a contact discontinuity. Of course, in view of the very low density of the disk gas (~ 10~^r3 "^^^g cm""^), the mass lost 

from the star in a single disk passage is bound to be very small, but so is the orbital energy; some TiRq/Mq ~ lO^r^^^ 

passages and ~ lO^rg^''^" yr elapse before circularization. A small minority of stars on orbits of low inclination and low 
eccentricity may encounter the disk at less than their own surface escape speed. Clearly the capture process presents further 
opportunities for interesting hydrodynamical, or radiation-hydro, simulations; Armitage et al. (1996) have simulated the 
passage of giant stars through a disk and find that much of the envelope can be stripped in a single encounter, but we are 
not aware of any hydrodynamic simulations for main-sequence stars. 

The second remark concerns the rate of accretion by the star once it settles into the disk. Syer et al. (1991) and 
Artymowicz et al. (1993) estimate this at the Bondi rate, 

M.B = ^^^^«5xlO-Vf/io yr-r^ 



11 



provided only that the Hih radius is smaher than the disk thickness, i.e. < 3OOr|M0. But accretion rate will be 
smaller if the gas is unable to cool sufficiently quickly; in other words, the rate should not exceed Eddington rate, 

^ = ^R^^ 1.0 X yr-i . 

Each of the first several doublings of the stellar mass is likely to require ~ 10^ yr, which may be gradual enough so that 
a gap will have formed by the time M* M/j. 

5. OBSERVABLE CONSEQUENCES 

Exotic and luminous as stars of ~ W^Mq would be, they could easily have escaped attention. Such stars radiate at their 

Eddington luminosity, but that is smaller than the Eddington limit of the quasar itself by the mass ratio M,/M ^ 10~^. 
Their effective temperature on the main sequence is « 7 x 10* K, based on the radius (eq. 8), which is determined 
essentially by virial considerations; a detailed atmospheric model might give a somewhat different photospheric radius 
and temperature. This is about twice the inferred Tes of bright QSOs (Malkan 1983). 

If the star has opened up a gap in the accretion disk, then this will affect the quasar's spectrum, particularly when the 
gap is in the innermost regions. In fact, because of the short viscous time of the inner disk, wc expect the gas inside the 
gap to be accreted relatively quickly, leaving an evacuated region between the central black hole and supermassive star, 
which is surrounded by what is effectively a circumbinary disk. The quasar will be relatively dim and red just before 
disruption or accretion of the star: indeed the source may not have been previously identified as a quasar, given the 
migration time of ~ 10^ yr. This is the time when any supermassive star would be most easily observed. For favorable 
inclinations, the star's spectral features would exhibit very fast orbital velocities of order several thousand kms^^ with 
periods greater than several days. The emission from the QSO disk, would gradually brighten and harden. Figure 3 
shows an approximate comparison between the spectrum of the supermassive star and QSO disk as the inner gas radius 
decreases. 

The merger of a supermassive star (or supermassive remnant) with the central black hole would be a strong source of 
gravitational waves. We estimate the "modified characteristic amplitude" for quadrupolc gravitational waves as defined by 
Finn & Thorne (2000) from the inspiral of a compact source with a mass equal to the isolation mass formed at lO^Rs from 
the central black hole (Figure 4). The trajectories in the frequency-amplitude diagram depend on the angular momentum 
of the central black hole, since the innermost stable orbit is much tighter for maximally rotating Kerr black holes: in 
the extreme case shown in Figure 4 with rotation parameter a=0.999, this radius is 1.18 GM/c^, while it is QGM/c^ for 
the Schwarzschild black hole. For the systems shown in Figure 4, a main sequence supermassive star would be disrupted 
before the differences due to angular momentum show themselves. Nevertheless, at and shortly after disruption, the 
waveform is bound to differ from that expected for an inspiraling black-hole companion. It would be interesting to try 
to predict these differences in some quantitative detail, as they might be the most direct way to confirm the existence of 
supermassive main-sequence stars in these disks. 

The dependence of the isolation mass on the mass of the black hole means that wc expect the supermassive stellar 
companions of more massive quasars to be more readily detected, although the signal is shifted towards lower frequencies. 
The cosmological redshift of the sources is important, as illustrated in Figure 4. The evolution of the gravitational wave 
signal from inspiral, should also be accompanied by a brightening and hardening of the quasar's electromagnetic radiation, 
as described above, and may be marked by an outburst if the star is disrupted. 

We can make a crude estimate of the upper limit to the event rate of supermassive star merger events from the known 
quasar population that may be seen by LISA. Using the quasar luminosity function of Boyle ct al. (2000), we estimate 
that the space density of "typical" quasars (within a couple of magnitudes of the knee of the luminosity function) is 
^ 10^^ Mpc^'^ in a flat Vlj^ = 0.7 cosmology. The comoving volume from 0.5 < z < 1. ■5 is 320 Gpc'' and from 1.5 < z < 2.5 
is 480 Gpc^ so there are 2.3, 3.5 x 10^ typical QSOs in these intervals, respectively, that can be seen over the whole sky. The 
supermassive star formation rate cannot be much greater than one per viscous time at rg.g. (and could be substantially 
less), which is about 1 x 10"'"' yr~^ per quasar. Thus the maximum event rate for LISA is of order a few supermassive 
star - black hole mergers per year. Of course the observed quasar luminosity function may underestimate the true number 
because of beaming effects. 

6. CONCLUSIONS 

We have argued that conditions in QSO accretion disks are likely to lead to the formation of massive stars. If the 
final masses of these stars are limited by the opening of a gap in the disk near the innermost radius where the disk is 
self-gravitating, then they would have ~ lO^M©. This is much more massive than any stellar object currently known. A 
key point is that the main sequence lifetime of these stars is much longer than their formation time, and modestly longer 
than the time they take to migrate to the black hole. Thus the end point of stellar evolution does not prevent the existence 
of these supermassive stars, and indeed the evolution may be interrupted only by tidal disruption and merger near the 
central black hole. A second important point is that once a star has formed in the disk with initial mass of order the 
Toomre mass, ^ IOOMq, then the subsequent accretion rates of surrounding disk material are much larger (~ Mq yr"""^) 
than the mass loss rates caused by thermonuclear instabilities (< 1O~'^M0 yr~^). 

Radiative feedback from the forming massive star (from luminosity generated either by nuclear burning or gravitational 
contraction) will have little influence while conditions in the disk and accretion flow are optically thick. Gas streamlines 
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will only be significantly perturbed as the stellar surroundings become optically thin, probably close to the time of gap 
formation, when the isolation mass has already been reached. 

Though they are very luminous, supermassive stars would easily be lost in the glare of the quasar. Periodic modulations 
of the light curve due to such stars would have periods ranging from a few years to a few days but with millimagnitude 
amplitudes. Even if the quasar was in a quiet state, perhaps because of the disruption of its disk by the supermassive star, 
then the light from the central regions of the host galaxy might dominate. The best candidates for periodicity searches 
would be quasars that appear to lack a blue bump, as this might be due to a gap opened in the disk by such a star. 
Space-based gravitational wave detectors are promising tools in the search for supermassive stars, particularly since the 
signal should also be accompanied by a rapid brightening of the quasar as the star is disrupted. 

Very recently, ideas similar to those presented here have hccii presented by Levin (2003). Although we have discussed 
our respective drafts, our two efforts have been independent. There is much we agree on, but there are some differences 
in assumptions and emphasis. Whereas Levin considers a disk that fragments into many stars of up to a few hundred 
solar masses, we consider the formation of a single dominant mass that accretes a significant fraction of the entire disk. 
Naturally we prefer our own scenario, but much more work will be required to decide which (if either) is correct. It is 
also possible that the outcome depends upon the initial conditions of the disk before fragmentation begins. Some insight 
into the early stages of fragmentation could be gained by numerical simulations of marginally self-gravitating disks; there 
have been many such studies, but none that we are aware of that include radiation-pressure dominance and radiative 
diffusion. Another difference between our studies is that and Levin considers parameters appropriate to the Galactic 
Center, whereas wc focus on bright and massive QSOs. Collin & Zahn (1999a) and Collin & Zahn (1999b) have also 
discussed star formation in AGN disks recently, but they too consider stars which, though massive by ordinary standards, 
are much less than the isolation mass considered here. 

Another area of theoretical uncertainty, which is perhaps even more challenging, is the nonlinear stability and evolution 
of extremely massive stars. In recent years, apart from some important contributions cited above, work on this topic has 
languished for lack of a clear astrophysical motivation. Convincing observational evidence for stars above ~ IO^Mq does 
not exist. But, as we have already stressed, the physical conditions in a quasar accretion disk are very much more extreme 
than those of a conventional star-forming region, and we have argued that the self-gravity and short timescales of these 
disks favor the rapid assembly of truly supermassive protostars. It would be very interesting to revisit the question of 
whether such objects can survive on the main sequence. 

We thank Yuri Levin and two anonymous referees for helpful discussions. JCT is supported by a Spitzer-Cotsen 
fellowship from Princeton University and by NASA grant NAG 5-10811. This work was supported in part by NSF grant 
AST-0307558 to JG. 
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Fig. 3. — Evolution of IO^Mq QSO accretion disk rest frame spectrum in arbitrary units (dotted lines) as the inner radius of optically thick 
emission migrates inwards from 500-Rs in steps of a factor of two. The outer radius is taken to be 1000i?s The spectrum (blackbody) from a 
main sequence, Eddington luminosity supermassive star of mass equal to the isolation mass is shown by the solid line. 
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Fig. 4. — Gravitational waves from the inspiral of a compact object of mass Miso that formed at r = 10^Rs,rB.g. (left, right panels), at 
redshift z = 1,2 (bottom, top panels; standard Cl\ = 0.7 cosmology assumed) from a supermassive black hole of mass IO'^Mq (lower lines) 
and 10* M0 (upper lines), as observed by LISA in a 1 year integration. For each case the "modified characteristic amplitude" for a harmonic 
of the waves (Finn & Thome 2000) from inspiral to a Schwaxzschild black hole is shown by the dashed line and to a Kerr bla<;k hole with 
rotation parameter of 0.999 by the dotted line. The square marks the point of disruption if the secondary is a main sequence star. The lower 
solid line shows LISA's rms instrumental noise level (Larson et al. 2000), averaged over the sky with normalization as in Finn & Thorne 
(2000) , while the upper solid line is the total noise including an estimate for the stochastic-background noise produced by white dwarf binaries 
(Bender & Hils 1997). 



